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ABSTRACT 

We present a VLT/UVES spectrum of a proximate sub-damped Lyman-a (sub-DLA) system at 
Zabs = 2.65618 toward the quasar Q0331-4505 (^q^o = 2.6785±0.0030). Absorption hnes of O I, Si II, 
Si III, Si IV, C II, C III, C IV, Fe II, Al II, and O VI are seen in the sub-DLA, which has a neutral 
hydrogen column density log A^hi = 19.82 ± 0.05. The absorber is at a velocity of 1 820 ± 250kms~^ 
from the quasar; however, its low metallicity [0/H]=— 1.64± 0.07, lack of partial coverage, lack of 
temporal variations between observations taken in 2003 and 2006, and non-detection of N V imply the 
absorber is not a genuine intrinsic system. By measuring the O VI column density and assuming equal 
metallicities in the neutral and ionized gas, we determine the column density of hot ionized hydrogen 
in this sub-DLA, and in two other sub-DLAs with O VI drawn from the literature. Coupling this 
with determinations of the typical amount of warm ionized hydrogen in sub-DLAs, we confirm that 
sub-DLAs are a more important metal reservoir than DLAs, in total comprising at least 6-22% of the 
metal budget at z « 2.5. 

Subject headings: intergalactic medium - quasars: absorption lines - galaxies: halos - galaxies: high- 
redshift 



1. INTRODUCTION 



Highly ionized gas is ubiquitous in galactic ha- 
los at high redshift. This is known from the de- 
tection of C IV and O VI absorption in damped 
Lyman-g (DLA, i.e. sys tems with logiVH j > 20.3; 



Ledoux et a l. 1998; W olfe fc Prochask a 2000,; iFox et all 
2007allbf) and sub-DLAs (systems w i th 19.0<logiy Hi < 



20 ^lDessauges-Zavadskv et al.l200"llPeroux et al.ll2003L 
l2007t l. These studies have found that C IV is present in 
every DLA where data exists, and that O VI absorption, 
tracing a hotter phase of gas, is present in at least 40% of 
DLAs. Detecti ons of C IV and O VI absorption have also 
been repor ted (jBergeron et al.l 119941 : iKirkman fc Tvtleil 
Il997lll999[) in Lyman Limit Systems (LLSs), a lower H I 
column density class of QSO absorber that may be trac- 
ing the extended outer regions of galactic halos. 

Physically, there are two principal processes that 
are expected to generate hot plasma in galactic ha- 
los. The first is star formation, which leads to 
supernovae and the subsequent production of hot, 
shock- h eated interstellar gas, possibly in the form of 
wind s (lOppenheimer fc Davel [20061 : iKawata fc RauchI 
120071: iFangano et al.l 120071) . The second is accretion, 
which results in the shock-heating of i nfalling intergalac- 
tic gas to tem peratur es of 10^~^K (Dave et al.l 1200 It 



iFang fc Brv"a3 2003; K ane et al.l |2005!). Given the low 
metallicities and low densities expected in these hot 
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halo environments, the cooling times will be very long 
(> lOiOyr for gas at T = lO^K, 0.01 solar metal- 
licity, and n = 10~''cm^'^), enabling metals to be- 
come locked up. Indeed, approximately half of the 
metals produced in stars by a redshift of two are 
yet to be observed (iBouche^et all I2005L l2006l l2007t 
ISommer-Larsen fc Fvnbcill2007f ). and hot halos represent 
a potential harbor for these metals. 

It has recently become apparent that sub-damped 
Lyman-a systems (sub-DLAs, with 19.0 < log A^hi< 20.3, 
also known as a Super-LLSs) play a sign ificant role in the 
metal budget (Peroux e t il]|2005l 120071 : iProchaska et all 
[2006; Kulkarni et al. 200^ but the contribution of the 
O Vl-phase in sub-DLAs has not been studied before; in- 
deed, we are aware of o nly one other publis hed detection 
of O VI in a sub-DLA (jSimcoe et al.ll2002t l. We present 
in this paper observations and ionization modelling of a 
sub-DLA showing a detection of O VI. 

2. OBSERVATIONS 

The spectra of QSO Q0331-4505 were taken with 
the Ultraviolet- Vi sual Echelle Spectrograph (UVES; 
iDekker et al.l l2000( l at the 8.2m Very Large Telescope 
Unit 2 (VLT/UT2) at Paranal, Chile during observ- 
ing runs in 2003 and 2006. The data reduction 
was performed using the UVES pipeline described in 
iBallester et al.l (|2000l ). The spectra were co-added in 
wavelength regions covered in both runs, and have a spec- 
tral resolution (FWHM) of 6.6 km s'^ (i?=45 000). 

A strong Lyman-a absorber exists in the spectrum at 
Zabs = 2.65618. This redshift is measured using the ve- 
locity of the strongest component of absorption in Si II 
A1304, and defines the velocity zero-point in the following 
discussion. Using a Voigt profile fit to the damping wings 
of the Lyman-a line, we find log A'hi = 19.82±0.05, iden- 
tifying the absorber as a sub-DLA. In Figure 1 we show 
the absorption line profiles of O I, C II, Si II, Fe II, Al II, 
Fe III, C III, Si III, Si IV, C IV, and O VI, which aU 
show detections in the sub-DLA. No C 11*, Al III, or N V 
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TABLE 1 
The Sub-DLA toward Q0331-4505 



Property 


Value 


Property 


Value 


^abs 


2.65618 




2.6785±0.0030=' 


d 


7.1Mpc'= 




1820±250kms-l 


log A^H I 


19.82±0.05 


log Ao I 


14.90±0.05 




76±2kms-l 




104±2kms-i 


log Afo VI 


14.43±0.02 


log Afc IV 


13.37±0.02 


logAfwarmHIl'* 


20.06±0.10 


[O/H] 


-1.64±0.07 


log AHotHII° 


>20.09 


A^HotHIl/A^HI 


>1.9 


log A^Total H II 


>20.38 


A'TotHIl/AfHI 


>3.6 



^ Zqao is an average of all the QSO emission lines measured. 

^ Estimated distance from the QSO, assuming Hubble flow (see text). 

Velocity width containing 90% of the integrated optical depth. 

-'^warm H II is derived from CLOUDY photoionization modelling to 
the low-ion column densities (see text). 

° The lower limits to A^Hot H II are found using the maximum allowed 
O VI ionization fraction (/ovi=0.2), and then calculating Ami = 
Afovi/(/oviO/H). 

absorption is seen. Both lines of the O VI doublet dis- 
play identical optical depth profiles over a velocity range 
of 175kms~^, verifying the genuine detection of O VI 
absorption, free from blending with the Lyman-a forest. 
Voigt profile fits were conducted using the VPFIT soft- 
ware package^, and are included in Figure 1. The fits 
to each ion were completely independent. A summary of 
the key properties of the sub-DLA is given in Table 1. 

We identify three principal components of low-ion ab- 
sorption, at —55, 0, and 20kms^^, which we name CI, 
C2, and C3, respectively. C2 is the strongest component. 
C3 is implied by the need to fit the positive- velocity wing 
of the low-ion lines, but its column density is relatively 
low, so we ignore it in the following discussion. CI and 
C2 are all seen in O I, Si II, Al II, Fe II, C II, C III, Si III, 
and Si IV, and C IV. We also identify two components of 
high-ion absorption seen in O VI and C IV, at —100 and 
— 60kms~^, which we designate HI and H2, respectively. 
The C IV seen near — 60kms~^ may have contributions 
from both H2 and CI. 

3. PROXIMITY TO THE QUASAR 

Proximate absorbers separated from the background 
quasar by dv <5 000 km s~^ [sometimes dv <3 000 km s~^ , 
where du=c|2;qso — Zabs|/(1 + ^abs)] are often removed 
from samples of intervening absorbers because they 
may be eit her ejected or ionized (or both) by th e 
QSO itself (|Wevmann et all HqtI iFoltz et all flQSl . 
We measured the observed wavelengths of the N V, 
O I, C IV, He II, and C III] QSO emission lines in 
the Q0331-4505 spectrum, and derived the QSO red- 
shift by accounting for the average velocities of these 
lines (relative to th e systemic redshift) as listed in 
iTvtler fc Fa^ (|1992f ). This yielded Zqso = 2.6785 ± 
0.0030^, implying the sub-DLA is separated from the 
quasar by a velocity of dw = 1 820±250kms~^. Assum- 
ing a WMAP 3-yea r cosmology (ffo=7 3kms~^ Mpc~^, 
r2M=0.24, » A=0.76:lSperge et al.ll2007D . we use the rela- 
tion given in lPhillipps et al! (|200S ) to find that the Hub- 
ble parameter at z = 2.65 is _ff(2.65)=257kms~^ Mpc~^. 
Therefore, assuming the sub-DLA motion is purely due 



^ Available at |http: / /www. ast.cam.ac.uk/ '^rfc/vpfit.html| 
® Improving the estimate of 2qsi-<=2.6 bv lMaza et al.l I I1993I) . 



to Hubble fiow (i.e., ignoring peculiar velocities), we find 
a separation of w 7.1 Mpc from the quasar. 

Despite the proximity to the quasar, several lines of 
evidence suggest this sub-DLA is not an intrinsic system 
arising near the central engine of the AGN: (1) the mctal- 
licity is too low (jPetitjean. Ranch, fc Carswelll994): (2) 
N V, which is expected in intrinsic systems ( HamannI 
|1997() , is not detected; (3) no excited states (e.g. C 11* 
or Si 11* ) are seen in ab sorption, as would be detect ed in 
intrinsic absorbers (e.g. ISrianand fc Petitieanl l2000): (4) 
there is no evidence of partial coverage of the background 
sourc e, which is found in intrinsic absorbers (jArav et al.l 
[2005h : (5) there is no evidence for temporal variability 
between our two observation epochs in 2003 and 2006. 

4. METALLICITY AND PHOTOIONIZATION MODELLING 

We use the O I/H I ratio to determine the metal- 
licity of the gas. Under interstellar conditions, O I 
is thought to f ollow H I closely due to charge ex- 
change reactions (|Field fc Steigma"iilll971h . implying that 
the O I/H I ratio accurately determines O/H with- 
out having to apply ionization corrections. Using a 
Voigt profile fit to O I A1302 and A1039, we find 
logA^(0 I) = 14.84±0.05 (total over the three compo- 
nents). Coupling this with logiV(H I) = 19.82±0.05, we 
find [O I/H I] = [O/H]=-1.64±0.07^ 

To determine the ionization balance in the photoion- 
ized gas, w e ran a grid of CL OUDY (version C06.02, 
described in lFerland et al.lfl998[ ) photoionization models 
to component C2 at different ionization parameter U = 
n^/nn. The models assume the gas exists in a plane- 
parallel, uniform slab, and take the metallicity [O/H] and 
A'hi as fixed inputs. We assume the fraction of the total 
H I column that exists in C2 is the same as the fraction 
of the total O I column in C2, which is 86%. To spec- 
ify the radiation field, we take an unattenuated quasar 
(QSO) spectrum, i.e. a power-law with the form oc 
at A < 1216 A (jHaardt fc Madaulll996l ) normalized 
by the measured QSO F-band magnitude of 17.9 and 
Zqso = 2.6785, which yield a total QSO luminosity at the 
Lyman Limit Lqi2=\.2 x 10^^ ergs~^ H z~^ when extrap- 
olatin g the F-band fiux as described in iGuimaraes et al.l 
(|2007( ). Making use of the QSO-absorber luminosity dis- 
tance of 7.1 Mpc, we then find that the flux of ionizing 
radiation from the QSO at the absorber is Fgi2=2.0x 
10^^^ ergcm"^ s""'^ Hz^^, lower than the extragalactic 
background (EGB) at z =2.65, which has f 912 =9.6 x 
10-21 erg cm-^s-^Hz-i (iHaardt fc Madaul[l99l . The 
QSO radiation field contains more high-energy photons 
than the EGB, boosting the production of Si IV and 
C IV, so the radiation field we use in the CLOUDY mod- 
elling includes both the EGB and the normalized QSO 
contribution. 

We found the value of \ogU that matches the 
Si IV/Si II column density ratio measured in C2. This 
value, logC/=— 3.32±0.10, also matches the observed O I 
column density. We then found that [C/0]=— 0.3 was 
necessary to reproduce the column densities of C II and 
(within a 2a error) C IV. Finally we checked that the 
model is consistent with the C III and Si III lower limits. 
This model is shown in Figure 2 and is characterized by 

^ Throughout this p aper we use th e solar oxygen abundance 
(0/H)q = 10-3-34 from lAsplund et al.l ll200i) . 
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Fig. 1. — VLT/UVES spectra of neutral, low-ion, and high-ion absorption in the sub-DLA at z^bs = 2.65618 toward Q0331-4505. This 
redshift, where the strongest component of low-ion absorption falls, defines the zero-point of the velocity scale. Our overall Voigt component 
fits to each profile are shown with solid red lines. Vertical dashed lines illustrate the positions of the three low-ion components (CI, C2, 
and C3) and two high-ion components (HI and H2). Note that in some panels, the profiles are shown over a limited y-range, for clarity. 
We have shaded the C IV and O VI panels to highlight the different structure in the high-ion absorption profiles. 



^ I I I I I ^ I I I I 
'-G2 Photoioniza^^on Model 

-In: [0/H] = -1.64; log 1=19.76 

^Out: log N^_^=20.00; x^"_= 0.64 [C/0] = -0.3 [Si/0]= 0.0 




log U 

Fig. 2. — CLOUDY photoionization modelling of component C2, 
assuming the radiation field has contributions from both the quasar 
and the extragalactic background. Each series of connected sym- 
bols shows the predicted column density of a given ion as a func- 
tion of ionization parameter (or equivalently, gas density, which 
increases to the left). The observations are shown in red, plotted 
at the density of the best-fit model. 

a hydrogen ionization fraction a;H=0.64, and an ionized 
hydrogen column density log7VHii=20.00. We assume 
the same conditions apply in CI and C3, giving a total 

log TVwarmH 11 =20.06. 

The ionization level we derive in the low-ion phase 
(64%) is no larger than the levels derived in inter- 
vening sub-DLAs by other authors (jProchaskal Il999t 
iPeroux et al.ll2007l ). reinforcing the idea that this sub- 
DLA does not show unusual properties because of its 
proximity to the quasar. The idea that many proximate 
DLAs (and by extension sub-DLAs) may be representa- 
tive of t he intervening popu lation has b een proposed be- 
fore (e .g. lMoller et al.lll998l ). though see lProchaska et all 
(1200% . 

5. A HIGHLY IONIZED PLASMA PHASE 



The O VI and C IV profiles are clearly different from 
the low-ion profiles, with the strongest high-ion compo- 
nent (HI) separated by 100 km s^^ from the strongest 
low- ion component. O VI and C IV show similar pro- 
files at —150 to — lOOkms^^, but the O VI absorption is 
smoother and broader than the C IV in the range —100 to 
km s^^ . This can be seen in Figure 3, where we compare 
the apparent column density profiles of the two ions. The 
difference in the profiles implies that the O VI and C IV 
are not fully co-spatial, i.e. the two ions trace different 
temperature regions of a multi-phase structure*. We de- 
rive a lower limit on A^h ii in the O Vl-bearing gas using 
A^Hot Hii = A^o vi/(/oviO/H), with an O VI ionization 
fraction /ovi — 0.2 (the maximum allo wed by coUisional 
ioniz ation models; see the Appendix in 'Trip p fc Sava"g^ 
I2OOOI ) and assuming that 0/H in the high-ionization gas 
is the same as 0/H in the low- ionization gas. This 
method yields logA^HotHii >20.09 (total over HI and 
H2). 

To explore how the O VI properties of our sub- 
DLA relate to other systems, we investigated the only 
intervening sub-DLA with a published O VI detec- 
tion: the syst em at Zabs = 2.439 toward Q1442-I-2931 
(jSimcoe et al.l (2002), and the previously unpublished 
O VI profile in the sub-DLA at Zabs=2.139 toward Q1037- 
270 (jSrianand fc PetitieanI [20011. The three sub-DLAs 
are described in Table 2. Following the same calculation 
as used here, the ratio A^HotHii/Ani takes the value >12 
in the Q1442+2931 sub-DLA, and >1.0 in the Q1037-270 
sub-DLA, vs >1.9 for the sub-DLA in this paper. From 
this we conclude that (i) in sub-DLAs, there are more 

^ Similar differences between the profiles of C IV and O VI 
have been found before in man y other environments, includ- 
ing the Milky Way halo JFox et a l. 2003) , the Large Magellanic 
Clou d IILehner fc HowM 120071') . LLSs (jKirkman fc Tvtleil [T997I . 
119991). D LAs ( Fox et alj [2007a|) , and some intergalactic absorbers 
HSimcoe et al.n200a V 
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Fig. 3. — Comparison of apparent column density profiles of 
C IV and O VI in the sub-DLA. Both ions show similar profiles at 
— 150 to — lOOkms"'^, but the O VI absorption is smoother and 
broader than the C IV in the range —100 to Okms~^. 

baryons in the O VI phase than there are in the form of 
H I atoms, and (ii) there is no evidence that radiation 
from the quasar is creating more O VI in the Q033 1-4505 
sub-DLA than in intervening sub-DLAs. 

6. CONTRIBUTION OF SUB-DLAS TO U 

The latest determination of the cosmic density in 
H I atoms in sub -DLAs is flfJ'i^DLA^j^ g ^ ;,^q-4 
([O'Meara at al.l I2007D . In order to fully account for 
the contribution of ionized gas in sub-DLAs, it is nec- 
essary to determine the average ionized-to-neutral ra- 
tio (A^TotalHIl/A'^Hl), where A'^TotalHII = ATwarmHII + 

A'^HotHii- The value of i Vwarm Hii/A^Hi has been inves- 
tigated bv iPeroux et al.l ()2007| 1. who report a mean ion- 
ization fraction of 0.68 in 26 sub-DLAs modelled with 
CLOUDY, corresponding to (iVwarmHii/A'^Hi)==2.1. The 
value of the iVnotHii/A'Hi is not well constrained: as dis- 
cussed in §5, the derived values in the three existing sub- 
DLAs with O VI are >1.0, >1.9 and >12. Assuming that 
the mean iVnotHii/A^Hi in sub-DLAs lies somewhere be- 
tween these measurements, we find that (A'^TotaiHii/A^Hi) 



[^I^ = (1 + (iVTotalHIl/A^Hl))f^Si 

between 6.1 x IQ-* and 2.2 x lO'^. 

The total contribution of sub-DLAs to the metal 
budget then follows according to: ^FiirfHii = 
(Z/H)r2|J'}?^j'}. Assuming a typical sub-DLA mctaUicity 
of ([Z/H]}=-0.8 (iKulkarni et alllMil . then nf^^]l^^^ 

is between 1.8 x 10~^ and 6.7 x 10~^, i.e between 6% and 
22% of all the metals produced by star formation up to 
z « 2 (using a val ue for the total Clz from star forma- 
tion of 3.0 X 10~^: lF^rara et aLll2005( ). For comparison, 
^zm+Hii ~ 6 X 10"'^ (including both the neutral and 
high ly ionized phases; iProchaska et al.l 120051 : iFox et all 
l2007al) . This calculation shows that globally there are 
at least three times more metal atoms in sub-DLAs than 
DLAs, owing to both the higher mean metallicity and 
the higher mean ionization level in sub-DLAs. 

There are several reasons why this estimate of flz in 
sub-DLAs could rise. First, recent detections of sub- 
DLAs with super-solar metallicity f | Khare et al.l 120041 : 
IProchaska et al.ll2006t [Peroux et alll2006f) may imply the 
mean sub-DLA metallicity is higher than —0.8. Second, 
the typical temperature of the O Vl-bearing gas in sub- 
DLAs may be above 3x 10^ K, implyi ng the O VI ioniza- 
tion f raction will be lower than 0.2 ( Sutherland Sz Dopital 
Il993| l . Finally, if the O Vl-bearing gas was generated by 
the accretion and shock-heating of intergalactic gas, the 
ionized gas may have lower metallicity than the neutral 
gas (increasing A'^HotHii)- Each of these effects would 
serve to increase i^z above the 6-22% level discussed 
here. 
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6 Plasma in a Sub-DLA 



TABLE 2 

Three Sub-DLAs with O VI Detections 



QSO 






Type 


log A^H I 


[O/H] 


log A^o VI 


AfHotHii/Ani 


Ref. 


Q033 1-4505 


2.674 


2.656 


Proximate 


19.82±0.05 


-1.64±0.07 


14.43±0.02 


>1.9 


1 


Q1442+2931 


2.669 


2.439 


Intervening 


19.78±0.30 


-1.82±0.40 


15.00±0.01 


>12.0 


2 


Q1037-270 


2.201 


2.139 


Intervening 


19.70±0.05 


-0.31±0.05 


Ril5.35±0.30 


>1.0 


3 
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